Hadronic Resonances from Lattice QCD 
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Abstract. The determination of the pattern of hadronic resonances as predicted by Quantum Chro- 
modynamics requires the use of non-perturbative techniques. Lattice QCD has emerged as the dom- 
inant tool for such calculations, and has produced many QCD predictions which can be directly 
compared to experiment. The concepts underlying lattice QCD are outlined, methods for calculat- 
ing excited states are discussed, and results from an exploratory Nucleon and Delta baryon spectrum 
study are presented. 
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MOTIVATION 

Patterns observed in measured spectra have repeatedly inspired fundamental break- 
throughs in modeling reality. The periodic table of the chemical elements led to the 
construction of the atomic model, and atomic spectroscopy led to the development of 
quantum mechanics. Later, the categorization of subatomic particles into the geometric 
patterns of the 'Eightfold Way' spawned the creation of the quark model and ultimately, 
the advent of Quantum Chromodynamics (QCD). 

It is critical to the development of particle physics to determine if QCD gives rise to 
experimentally observed nuclear physics phenomena. The QCD spectrum is a particu- 
larly desirable calculation not only because it would provide a check of the theory, but 
also because it could lead to a much deeper understanding of the physics of the strong 
interaction. Additionally, because its underlying quark degrees of freedom are color con- 
fined, QCD's low-energy behavior must be inferred from the properties of its composite 
hadronic resonances. 

One of the main goals of the Lattice Hadron Physics Collaboration (LHPC) is to cal- 
culate the QCD hadron spectrum from first principles using the lattice QCD formalism. 
Such a calculation would shed light on many interesting phenomenological issues such 
as the nature of the Roper resonance, the missing quark model resonances, and the prop- 
erties of hybrid and exotic hadronic resonances. 
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LATTICE QCD AS A TOOL TO STUDY RESONANCES 



The need for the lattice formalism 

In QCD the effective strength of the strong nuclear force is controlled by the renor- 
malized coupling between the quarks and the gluons. This quantity is dependent on the 
energy scale of the interaction and becomes weak at high-energies, leading to the prop- 
erty of asymptotic freedom. This is consistent with the results of deep inelastic scat- 
tering experiments involving high momentum transfer interactions. At these scales, the 
application of perturbation theory has led to predictions which agree spectacularly with 
experiment. 

In contrast, the coupling becomes strong at lower energy scales. This is consistent 
with color confinement, and leads to a rich spectral structure. This strong coupling 
precludes any perturbative expansion (in the coupling) for low energy quantities such 
as the spectrum. Fortunately, non-perturbative calculations are possible using a space- 
time lattice regulator in conjunction with the functional integral framework [1]. The use 
of a lattice not only defines a measure for the functional integral, but also introduces an 
ultraviolet cutoff via the inverse lattice spacing aT x . In the following we survey several 
key concepts which enable lattice practitioners to calculate the QCD spectrum. 



The physics of lattice QCD 

Several expository books on lattice QCD are available [2]. Here we review those 
features of the formalism relevant to spectroscopy, using scalar fields for simplicity. In 
the continuum, quantum fields <j>(x;t) are defined over all spatial points x at each time t, 
and are related at different times through the time evolution operator. 

zxp[-M(t-t% (1) 

where H is the Hamiltonian of the theory. The lattice formulation uses fields restricted 
to the sites of a 4-dimensional hypercubic lattice. The lattice fields are denoted <j> T (x), 
where x is the spatial index now restricted to values on a 3-dimensional cubic lattice. The 
discrete index x does not represent the standard (Minkowski) time t, but instead denotes 
imaginary (Euclidean) time. The use of Euclidean time implies that the collection of 
fields on one time slice is related to the collection of fields on a neighboring time slice 
via the transfer matrix: 

T = exp[-aH], (2) 

where a is the lattice spacing between neighboring time slices. Using the exponential 
Euclidean transfer matrix allows us to utilize many computational techniques from 
statistical mechanics not available when using the oscillatory Minkowski time evolution 
operator. If we use a finite temporal extent with periodic boundary conditions, we have 
the system shown in Fig. 1. 



-aH 



-aH 



-aH 



e 



e 



e 



FIGURE 1. A graphical representation of the system lattice QCD simulates: N time slices of quantum 
field variables coupled by the transfer matrix T = exp(-aH). Periodic temporal boundary conditions 
identify the fields on the first (t = 0) time slice with those on the the last (t = Na) time slice. 



If we perform an integral over the possible values of the field variables on all of the 
time slices {0 T }, we get the partition function of the theory: 



where we have used the completeness relation / d0 T |0 T ) (0 T | = 1. 

QCD considerations. To extend this formalism to QCD, we introduce quark field 2 
variables on the sites of the lattice, and connect them with gauge link variables repre- 
senting the gluonic degrees of freedom. The gauge links provide a connection between 
adjacent sites that maintains color gauge-covariance. We use a finite temporal extent 
/3 = N r a, and render the number of integrals on each time slice finite by working in a 
finite periodic spatial volume L 3 = N s a , where N T and N s are the number of temporal 
and spatial lattice sites, respectively. The temporal extent /3 in Eq. 4 is analogous to the 
inverse temperature in statistical mechanics: 



and thus the lattice formalism can be used to calculate the QCD Equation of State 
and investigate such systems as the quark-gluon plasma [3]. To determine the QCD 
spectrum, we are interested in the zero-temperature limit of the partition function, and 
will work with large temporal extents j8 — > °°. It should be noted that Minkowski time 
does not appear anywhere in this formulation. This does not cause difficulties because 
a theory's spectrum is an equilibrium (i.e. time-independent) property. It is possible, 
however, to define lattice formulations which allow the investigation of non-equilibrium 
properties of QCD such as transport coefficients [4]. 
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j8 = N t a ~ 1 / (temperature) , 



(5) 



To maintain the definition of the trace, (fermionic) quark fields use anti-periodic temporal boundary 
conditions while (bosonic) gauge links use periodic temporal boundary conditions. 



The spectral representation of correlation functions 



To access the spectrum, we define a correlator C^\z) for < z < /3 by consider- 
ing correlation between time-ordered source and sink operators defined on time slices 
separated 3 by z (c.f. Fig. 1 and Eqs. 3-4): 



C0>(T) = <tf(T)tf(O)) = -Tr 
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We may express the traces in terms of the energy eigenstates: 
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Finally, we take /3 >> z to access the zero-temperature physics: 
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We may decompose the zero-temperature correlator C(t) defined in Eq. 9 into its 
spectral components by inserting a complete set of energy states \k): 
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The derivation of Eq. 12 neglects boundary conditions, and is therefore accurate only in 
the infinite- volume limit. In this limit we see that C(z) decays as a sum of exponentials 
with decay constants given by the energy levels accessible by the application of 6 to 
the vacuum state of the system. For large Euclidean time separations z, we see that 
the correlator is dominated by the energy gap between the vacuum and the first excited 
state. It is crucial to design operators having | (k\ ^|0) | 2 large for the states \k) of interest, 
and small for the other (contaminating) modes. An overview of the LHPC's operator 
construction method can be found elsewhere in these proceedings [5], and more detailed 
descriptions may be found in [6, 7]. 



IMPROVING SIGNAL QUALITY 

In practice, the correlator C(z) is estimated using the Monte Carlo method [8], and will 
consequently have an associated uncertainty for each value of z. It can be shown [9] that 



3 The temporal boundary conditions may be used to shift the source operator to time slice and the sink 
operator to time slice T. 
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FIGURE 2. Left: a local operator's effective mass signal. Quark field smearing drastically reduces the 
operator's coupling to short-wavelength contaminating modes. Right: an extended operator's effective 
mass signal (after quark field smearing). Gauge link smearing strongly attenuates the noise coming from 
the stochastic evaluations of the gauge links used in constructing the operator. 



the signal-to-noise ratio for baryon correlators decays exponentially with x. Recalling 
the expression for the correlator given in Eq. 12: 

oo 

C(T)= 1 £\(k\W\0)\ 2 e- E ^ 
k=i 

we see that if the operator & couples strongly to several states \k), the first excited energy 
level Ei may not have a chance to dominate before the signal is lost in the noise. 

We may define the effective mass function, which becomes the lowest energy (in 
lattice units) at large values of v. 

aM{%) =ln[C(r)/C(T + a)] ^ aE x . (13) 

Effective mass plots, such as those shown in Fig. 2, are a useful tool for visualizing and 
evaluating an operator's correlator signal according to two key criteria: (1) excited state 
contamination, and (2) noise. 

Fig. 2 demonstrates two powerful ways to improve signal quality: quark field smear- 
ing and gauge link smearing [10]. The smearing process replaces the variables at each 
location with a suitable form of a weighted local spatial average. It has been found 
that quark field smearing drastically reduces the coupling of operators to the short- 
wavelength contaminating modes of the theory, at the price of a modest increase in 
noise. The primary source of noise in an operator, especially in extended operators using 
quarks which are covariantly displaced from one another, is the presence of stochasti- 
cally updated gauge link variables. The application of link smearing strongly attenuates 
the noise of the operator. 



EXTRACTING EXCITED RESONANCES 



Using the variational method to extract excited states 

Instead of one operator &, we may use a basis of operators { ff a } to define a correlator 
matrix: 

(W^ b (0)) ^ (01^-^10) =Cab{T). (14) 

We may view the quantities C^(t) as matrix elements between states in the trial basis 
\ff a ) = a \O). Once we have estimates of these elements, we may apply the variational 
method [11] to define a new operator basis 0, = £ a & a v a U an d may choose the coeffi- 
cient vectors v,- to diagonalize e ; in the subspace spanned by our trial basis: 

v]c(r) Vj = (0|®,-[c- flT ]®y|0), (15) 
= Oif*W (16) 

This orthogonality implies that the diagonal matrix elements v?C(t)v/, the so-called 
principal correlators, will be asymptotically dominated by different energy levels, as 
shown in Fig. 3. One refinement of the method is to relax the trial basis by operating 
upon it with exp(— Htq) for some small relaxation time To- Thus, by diagonalizing 

C-WWCWC-Wfa), (17) 

it can be shown [12] that one is working (formally) with the basis of trial states: 

\^ a (T )) =cxp(-HT /2)W a \0). (18) 

By relaxing the trial basis, we obtain a better overlap with the subspace spanned by the 
low-lying states of interest. In practice, this method is complicated by the presence of 
noise in the estimates of C a },{x). 




FIGURE 3. A toy theory with three energy levels shown as dashed horizonal lines. Two operators are 
used to form a 2 x 2 correlator matrix, and the effective masses associated with the diagonal elements 
C aa (?) are shown on the left. Diagonalizing the correlator matrix yields the principal correlators v-C(t)v,-, 
whose effective masses are shown on the right. This method can be used to extract multiple excited states. 
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FIGURE 4. Left: The low-lying Nucleon (top) and Delta (bottom) baryon spectra as determined by 
200 quenched ("zero flavors of dynamical quarks") configurations on a 12 3 x 48 anisotropic lattice 
with a s ~ 0.1 fm and a s ja x ~ 3.0. The quenched approximation neglects the effects of dynamical sea 
quarks, and corresponds to neglecting quark loops in diagrammatic calculations. Our pion mass for this 
exploratory study was approximately 700 MeV. The vertical height of each box indicates the statistical 
uncertainty in that estimate. Right: The corresponding spectra as determined by experiment [13] and 
projected onto the space of lattice spin-parity states. While there are infinitely many continuum J p 
quantum numbers, there are only a finite number of irreducible representations of the spinorial lattice 
rotation group O® : G\ g , H g , and G2 g having even parity, and G\ u , H u , Gi u having odd parity. Consequently, 
continuum P states will appear on the lattice as degenerate levels in one or more lattice spin-parity 
channels. For example, the continuum Dis(1675) J p — | + resonance appears on the lattice as a pair of 
degenerate levels appearing in the H u and G21, spin-parity channels. 



RESULTS AND OUTLOOK 



Lattice QCD calculations of the Nucleon and Delta baryon spectra [5, 7] are presented 
in Fig. 4. Comparison is made with experimental results [13] by subducing continuum 
J p quantum numbers onto their discrete lattice spin-parity counterparts [6] (even parity: 
Gi g , H g , Gi g , odd parity: G\ u , H u , G2J. Although this exploratory study is quenched 
(neglecting the effects of dynamical sea quarks) and performed at an unphysically high 
pion mass of approximately 700 MeV, several interesting patterns are seen. 

The quenched spectra reproduce the isolated even-parity states corresponding to the 
proton and the A ++ (1232). We also see a band of low-lying odd-parity states in each 
spectrum matching the corresponding low-lying experimental odd-parity bands. In con- 
trast to experiment, we find that the first excited even-parity states corresponding to the 
Roper and the P^{\600) resonances lie above the first odd-parity band of states. Also, 
the splitting between the first excited even-parity band and the lowest-lying odd-parity 
band of states is significantly larger in the quenched spectrum. These discrepancies may 
be due to quenching, sensitivity to the quark mass, finite-volume effects, discretization 
artifacts, or poor interpolation by our three-quark operators. Unquenched runs currently 
underway at multiple volumes and pion masses will attempt to resolve these issues. 

In conclusion, the lattice formulation is a powerful method for calculating non- 
perturbative quantities used in fundamental tests of QCD. The exploratory results pre- 
sented here demonstrate clear progress toward the long-standing goal of determining the 
QCD resonance spectrum from first principles. Ongoing studies using improved opera- 
tors, variational methods, and dynamical sea quarks will continue to reveal the patterns 
present in the QCD spectrum. 
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